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ABSTRACT 

We use a modified outer gap model to study the multi-frequency phase- 
resolved spectra of the Crab pulsar. The emissions from both poles contribute 
to the light curve and the phase-resolved spectra. Using the synchrotron self- 
Compton mechanism and by considering the incomplete conversion of curvature 
photons into secondary pairs, the observed phase-averaged spectrum from 100 
eV - 10 GeV can be explained very well. The predicted phase-resolved spectra 
can match the observed data reasonably well, too. We find that the emission 
from the north pole mainly contributes to Leading Wing 1. The emissions in 
the remaining phases are mainly dominated by the south pole. The widening of 
the azimuthal extension of the outer gap explains Trailing Wing 2. The com- 
plicated phase-resolved spectra for the phases between the two peaks, namely 
Trailing Wing 1, Bridge and Leading Wing 2, strongly suggest that there are 
at least two well-separated emission regions with multiple emission mechanisms, 
i.e. synchrotron radiation, inverse Compton scattering and curvature radiation. 
Our best fit results indicate that there may exist some asymmetry between the 
south and the north poles. Our model predictions can be examined by GLAST. 
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Introduction 



It is generally believed that the phase-resolved spectra can provide the most detailed 
information about the structure of the pulsar magnetosphere, accelerat ion mechanism , and, 
pair creation and rad i ation processes in the outer magnetosphere. iRomanil (119961 ) and 



Daughertv fc Hardind f|1996T) firstly calculated the phase-resolved sp e ctra o f the Vela pul 



sar. Then, ICheng. Ruderman fc Zhang! (120001 ) and IZhang fc Cheng! (120021 ) calculated the 
phase-resolved spectra of the Crab pulsar in the 7-ray and the X-ray regimes separately. It 
will be interesting to find a general scenario that can produce both the 7-ray and the X-ray 
regimes in the phase- resolved spectra of th e Crab pulsar. Our preliminary results have been 
reported in the 2006 COSPAR conference (IJia et al.l 120071 ). We have continued to study the 
phase-resolved spectra of the Crab pulsar in mo re details. Here, we present our recent results 
with a different parametric fitting (ITangi 120071 ). 



According to ICheng. Ho &: Ruderman! (Il986al : Il986bl : hereafter CHR), the outer gap 
starts at the null charge surface, ends at the light cylinder, is bounded below by the last 
closed field line and is bounded above by a layer of electric current which replenishes charges 
to the open field lines outside the gap to maint ain a steady charge density, the Goldreich- 
Julian charge density (jGoldreich fc Julian! 1 19691 ) . 
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The charge depletion within the outer gap due to global flows of charged particles causes a 
large electric field along the magnetic field lines so that inside the gap. This gap 

thus acts as an accelerator to boost the charged particles to relativistic speeds. Through a 
cascade process, high energy 7-ray phot ons and e^ pairs are produced . Rece ntly, the classical 
outer gap model is being challenged by lHirotani. Harding fc Shibatal (120031 ). By solving the 
set of Maxwell and Boltzmann equations, they find that a current at nearly Goldreich- Julian 
rate can shift the position of the inner boundary of the outer gap. Therefore, we adopt a 
modified version of the CHR outer gap model such that the inner boundary of the outer gap 
is shifted inwards. 

The photon emission mechanism starts with curvature radiation of the accelerated 
charged particles in the gap. The emission direction is tangent to the local magnetic field 
lines. As a photon escapes, it may encounter a low energy photon. The low energy photon 
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may be a thermal photon from the stellar surface or a magnetospheric soft synchrotron pho- 
ton emitted by the secondary e pairs which are created by the curvature photons from the 
inner field lines. The primary curvature photons will then be converted into the secondary 



e 1 * 1 pairs via photon-photon pair production, i.e. 7 + 7 



+ e or 7 + X 



e+ + e~ 

inside and outside the gap. As pointed out in lCheng. Ruderman fc Zhang! (120001 ) . although 
pair production inside an outer gap is limited to a small region, pair production outside 
the outer gap can cover a much wider area because the synchrotron photons produced by 
the secondary pairs are more abundant than the thermal photons from the stellar surface. 
Although the secondary synchrotron photons cannot get into the outer gap due to the field 
line curvature, they can convert most of the primary curvature photons from the outer gap 
into secondary pairs. Sin ce the synchrotron radiatio n is beamed to a small angle by the rel- 
ativistic beaming effect (IRybicki &: Lightmanlll979l ). synchrotron photons will also be seen 
as more or less tangent to the field lines if observed. 

Besides curvature radiation and synchrotron radiation, inverse Compton scattering is 
another important radiation mechanism in the neutron star magnetosphere. It occurs when 
a fast moving electron or positron collides with a photon and net energy is transferred from 
the particle to the photon. In the far region of the magnetosphere, the relativistic particles 
collide with the soft synchrotron photons through the inverse Compton scattering process. 



Chiang fc Romanil (119921 . Il994l ) and iRomani &: Yadigaroglul (119951 ) calculated the light 



curves by considering a single outer gap with the photon emission to be beamed to the outside 
alone in order to avo id the multiple-peak feature that do es not occur in true observation data 
of the Crab pulsar. ICheng. Ruderman fc Zhang! (120001 ) pointed out the lack of a reason to 
explain the ignorance of the incoming photons (the photons that are beamed towards the 



star). However, in Section [27T 
the outward emission (see also 



we estimated that inward emission is much fainter than 



Cheng. Ruderman fc Zh ang 20 00|). Furthermore, since the 



previous models by IRomani fc Yadigaroglul (119951 ) and ICheng. Ruderman fc Zhang! (120001 ) 
considered emission beyond the null charge surface only, they suggest that the observer 
can measure the photons from one gap only, that is, a single-pole model is considered in 
the canonical model. However, if the gap extends below null charge surface, the photons 
originating from another gap become measurabl e by the obs erver as well. In this paper, 



theref ore, although we follow IChiang Romanil (119921 . 119941 ) and IRomani fc Yadigaroglu 
(119951 ) to calculate the light curve, we consider both gaps, that is, a two-pole model is 
examined. 

In this paper, we calculate the pulse profile and the phase-resolved spectra for the Crab 
pulsar with an outer gap accelerator model. Because the Crab pulsar is one of the brightest 
7-ray sources in the sky, the detailed observation for the pulse profile and the phase-resolved 
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spectra have been obtained. These provide useful information to the study of the non- 
thermal processes in the pulsar magnetosphere. Since the Crab pulsar is still young and is 
believed to have a thin gap, calculating the contribution from one layer suffices. By fixing 
an inclination angle and an observer angle, the photon emission locations that produce the 
light curve are used to calculate the phase-resolved photon spectra which are compared with 
the observation data. For the Crab pulsar, we adopt R = 10 km for the stellar radius and 
B p = 3.8 x 10 12 Gauss for the stellar magnetic field strength. 

Another major feature in this paper is the relaxing of the assumption that all curvature 
radiation has been converted into e pairs as in the previous studies. In fact, by letting trace 
amount of the curvature radiation, which is emitted far away from the star, to survive and 
escape, we can obtain better fitting curves for the phase-resolved spectra of the Crab pulsar, 
especially in the phases of Trailing Wing 1, Bridge and Leading Wing 2. 

In Section [21 we review briefly the emission geometry in the magnetosphere. Then we 
discuss the location of the inner boundary of the outer gap and produce the light curve for 
the Crab pulsar. In Section [31 we discuss the electric field component along the magnetic 
field lines and the three major emission mechanisms, namely synchrotron radiation, inverse 
Compton scattering and curvature radiation. We argue that although most of the primary 
curvature photons have been converted to secondary pairs to produce synchrotron photons, 
some curvature photons which are emitted far away from the star can survive the cascade 
process and escape. They contribute to the peak in the high energy regime at several GeV 
of the phases Trailing Wing 1, Bridge and Leading Wing 2. In the final section, we conclude 
our results and discuss briefly about the justification of our assumptions. 



2. Theoretical light curves of neutron stars with Crab parameters 

2.1. Emission geometry 



To calculate the light curves and the spectra, we adopt the rotating dip ole field in the 



magn etosphere. For a rotating dipole, the local magnetic field B(r) is given by (jCheng. Ruderman &: Zhang 
2000) 



B 



3 2 7" 



A* M A* 



c 2 r 



where /x = /i (sin a cos fltx + sin a sin Qty + cos az) is the magnetic moment vector, f is the 
radial unit vector and a is the inclination angle. 



We calculate the polar cap edge by using (x , y , z ) = (R p cos(j) p , R p smcj) p , sjR 2 — R p ) 
as the initial trial values for the computer program. R p = R\J R/ Rl is the polar cap 
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radius of an aligned static dipole where R is the stellar radius and Rl = c/Q is the radius 
of the light cylinder. <p p is the azimuthal angle about the magnetic axis. We call it the 
polar cap angle. Then we employ the Runge-Kutta method to trace out the field lines. 
In subsequent iterations, we try to find the scaling factors, a , such that (x' , y' , z' ) = 
(a x , a y , \ R 2 — a lR 2 ) correspond to the footprints of the last closed field lines, i.e. the 



boundary of the polar cap. In Figure 2 of ICheng. Ruderman fc Zhang) (120001 ) . they have 
shown that a is p -dependent. The polar cap of a rotating dipole is not circular in shape 
especially for those having large inclination angles. A three-dimensional view of the last 
closed field lines is shown in Figure HJ Next we define another scaling factor a± to denote the 
footprints of the open field lines as (x,y,z) = (aix' ,aiy' , a/ R 2 — (x 2 + y 2 )). Here a\ — 
represents the magnetic pole and a x = 1 represents the last closed field lines. 



Since the star is rotating, aberration occurs along the line of sight. With (3 
we have 



\r x Q\/c, 



u, 



(u^ + (3c) 
(1 + (3u^/c) 



P 2 



[1 + f3u4,c) 
[1 + Pu^c) 



(3) 



where u% and i4 for i = r, 9 and <fi are the emission direction in the co-rotating and observer 
frames, respectively. Choosing the rotational a xis as the z-axis, in the observer frame, the 
polar angle from the rotational axis is given by (lYadigaroglulll997l ) 



it, 



cosC = -f. 
u' 



(4) 



where ( is the viewing angle with ( = 0° when the star is viewed directly above its rotational 
axis and ( = 90° when the star is viewed over the stellar equator. 

Comparing to the photons emitted at the centre of the star, a photon emitted at any 
particular location r will take less time to travel to the light cylinder. The phase difference 
due to the travel time is given by A$ = — r ■ v! /Rl- Therefore, we have the phase angle $ 
in the observer frame (j Yadigaroglul 1 1 9971 ) to be given by 



$ = -<//- 



r ■ u 



Ri 



(5) 



where —<$' = — cos 1 {u' x /u' xy ) is the azimuthal angle in the observer frame. Choosing f2-/x 
plane to be the x-z plane, u' xy is the length of the projection of v! on the x-y plane. 
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We would li ke to remark that we ass u me th e dipole field in the corotating frame. On 
the other hand, iTakata. Chang Cheng (120071 ) assumed the rotating dipole field in the 
observer frame. With the difference in the magnetic field configuration in the observer 
frame, the emission direction near the l ight cylinder is azimutha l direct ion for the present 
case and nearly radial for the case in ITakata. Chang fc Chengl (120071 ) . Essentially, this 
difference appears because the poloidal magnetic field dominates the toroidal field near the 
light cylinder in the observer frame for the pres ent case, and the poloida l and t he toroidal 
fields are comparable to each other for the case in ITakata. Chang fc Chengl (120071 ) . Although 
there is a large differen ce in emission direction near th e ligh cylinder, the pulse profiles do 
not change very much (ITakata. Chang &: Chengl 120071 ) because the radiation very close to 
the light cylinder mainly contribute to the bridge phase. As we will see later, however, the 
calculated ph ase-resolved spectra in the pres ent case explain the observations better than 
the results in ITakata. Chang fc Cheng (120071 ). Therefore, the present results suggest that 
the radiation diretion near the light cylinder seem to be in the azimuthal direction rather 
than in the radial direction. Recen t particle simulati o n for the global structure of the charge 
separated magnetosphere done by IWada fc Shibatal (120071 ) also indicated such behavior of 
the emission direction. On the other ha nd, the time- dependent force free relativistic MHD 
solution obtained by (jSpitkovskyl 120061 ) indicated the radial motion of the particles near 
the light cylinder. Furthermore Bucciantini et al. (2006) have considered a more general 
relativistic MHD approach for rotating pulsars a nd their solutions asymptotically approach 
the force-free ones similar to that obtained by (jSpitkovskyl 120061 ) in the high magnetized 
wind case, which is closed to our case. Therefore the global structure of the magnetic field 
is still an open question. 

In this paper, we neglect the contribution from the inward emissions with the following 
reasons. In general, since the charged particles accelerate only within the gap and they lose 
energy during the cascade process, the incoming charged particles and hence the incoming 
photons cannot have an energy exceeding r j e m e c 2 which is the energy of a charged particle 
immediately after it leaves the gap. 7 e ~ 10 7 is the local Lorentz factor of the charged 
particle upon leaving the outer gap (cf. Equation [8]). On the other hand, charged particles 
within the gap are being accelerated con tinuously and will gain an e nergy of eV gap with 



V gap « 6.6 x 10 12 f$B l2 P- 2 V ~ 10 15 V (iRomani fc YadigaroghJ Il995h where f « 0.2 is 



the average value of the local gap size at Rl- Therefore, we can roughly estimate the ratio 
of the intensity of the radiation due to the incoming photons to the radiation due to the 
outgoing photons by (N gap ^ e m e c 2 ) / (Ng ap eVg a p) ~ 0.5%. As a result, when we compute the 
light curve, we neglect the contribution from the incoming photons. 



By considering the radiation to be emitted tangent to the magnetic field lines in the 
co-rotating frame, we project the radiating points onto the £ — $ plane. Figure [2] shows 
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the photon emission pattern for the inclination angle a = 50° and a\ = 0.97. Here we 
assume the emission region to extend from the stellar surface to the light cylinder and 
assume a symmetry between the north and the south poles. In other words, when there is 
a photon emitting with (£, $) from the north pole, there is another photon emitting with 
(180° — C, 180° + $) from the south pole as well. In Figure [2j the grey lines correspond to the 
outgoing photons emitted from the north pole, the pole which is making an acute angle with 
the rotational axis. The black lines correspond to the outgoing photons emitted from the 
south pole. For example, for an observer at a viewing angle smaller than 90°, the emission 
region by the north pole (grey) corresponds to the radiation emitted in the region between 
the inner boundary (stellar surface) and the null charge surface, and the emission region by 
the south pole (black) corresponds to the radiation emitted in the region beyond the null 
charge surface. 

The viewing angle which is the angle between the observer and the rotational axis can 
be set at a certain value Co- Then we can measure the number of photons traveling in this (o 
direction and produce a theoretical light curve. However, before we move on to produce the 
theoretical light curve, we need to mention that the radiation is, in fact, emitted within a 
finite emission cone of half- angle <p{r) instead of simply tangent to the field lines. The criteria 
for counting a photon becomes ( — (p(r) < (o < ( + <p(r). This effect can be understood 
in a geometrical point of view. Figure [3] shows a close-up illustration of two field lines that 
approximate two concentric circles. h(r) is the local thickness of the outer gap, s(r) is 
the local radius of curvature of the field line and A(r) is the pair creation mean free path. 
According to Figure [31 since the secondary pairs are produced just above the outer boundary 
of the outer gap accelerator, we can estimate the pitch angle of the new born pairs as 

. 2 2f(r)R L 

8111 ^ (r) = s(r) ' ( } 

where f(r) is the fractional gap thickne ss defined by f(r) = h{r)/Ri,. The self-sustained 



outer gap model of lZhang fc Cheng! (119971 ) estimated the fractional gap thickness as f(Ri/2) 



5. 5 p26/2l£-4/7 which 

is ~ 0.04 for the Crab pulsar and ~ 0.13 for the Vela pulsar. By 
considering the conservation of magnetic flux along a field line, with 1-D approximation, 

3 



f(r) = f(Ri) Cme) 2 f Cheng. Ruderman fc Zhang 2000 ). Since s(r) = \frRl in static 



dipole approximation, sin<£>(r) = sin(p(Ri) 



1/2 



2.2. Inner and outer boundaries of the outer gap 



As the electro dynamical studies have shown (Hirotani. Harding &: Shibatall2003 ; Takata. Shibata &: Hi: 



2004 ; iTakata. Chang &: Chena 120071 ; iHirotanil l2006bl ) . the inner boundary of the outer gap 
accelerator is shifted inwards from the null charge surface with an increase in the current 
through the gap. In fact, if there is no current injection from the inner and the outer bound- 
aries, the inner boundary will be located at a position on which B z /B = j g is satisfied, 
where j g is the current density in unit of QB /2tt carri ed by the pairs created in the ga p 
and is constant along the field line for the steady state (ITakata. Shibata fc Hirotanill2004l ). 
For example, if no current is created in the gap (j g = 0), the inner boundary is located 
at the position where B z /B = is satisfied, that is, on the null charge surface. On the 
other hand, if j g ~ cos a on a particular magnetic field line, the inner boundary on the field 
line is located at the stellar surface, where B z /B ~ cos a is satisfied. We expect that the 
created current de nsity is proportional to the pair-cre ation rate, which depends on the radial 
distance as r -3 / 8 ( Cheng. Ruderman &; Zhangl l2000i ) . Since most of the pairs are created 
around the null charge surface in the outer gap, we may be able to relate the created current 
density with the radial distance to the null charge surface of each last closed field line as 
jg(4>p) — jg(®)l r n(0) /r n (4>p)]~ 3 ^ 8 , where j g (0) and r n (0) are, respectively, the created current 
density on and the radial distance to the null charge surface on the last closed field line 
with the polar cap angle <p p = 0°. We know the location of the inner boundary with the 
azimuthal angle if we would estimate the created current density j g (0). As demonstrated 
by the electrodynamical studies, however, the current structure in the gap is very sensitive 
to the gap geometry such as trans-field thickness and the longitudinal width. Since the 
gap geometry in the p ulsar magnetosphere should be determined by the global condition 
( IWada fc Shibata! 120071 ) and since there is no study for the 3-dimensional magnetosphere of 
an inclined rotator, we deal with the created current j g (0) by using some model parameters. 
Figure H] summarizes the variation of the radial distances to the inner boundary and the 
null charge surface on the last closed field lines against the polar cap angle of the field lines 
around the magnetic axis. From Figure HJ we can see that only a small current is created 
around P ~ 180°, so the outer gap must be less active there. In this paper, we constrain 
the width of the polar cap angle with the field lines on which the pair creation mean-free 
path, A(r) ~ [2s(r)f(r)R L 



1/2 



2/ 1 / 2 (i?i/2)r, at the null charge surface of the active field 



lines is estimated to be shorter than R^. This condition produces an azimuthal extension of 
the outer gap of A<p p ~ 250°. 



as 



For the outer bounda ry, the position should be determined with the global model such 



Wada fc Shibata! (120071 ) with current. For the present local model, the position is a free 



parameter and we put it at the light cylinder because we assume that the emissivity of the 
curvature radiation is declined rapidly near/beyond the light cylinder and/or the radiations 
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beyond the light cylinder are beamed out of line of sight due to the magnetic bending. We 
may put the outer boundary inside of the light cylinder. However, the resultant pulse profile 
and the spectra does not change very much unless the outer boundary is located close to 
the null charge surface, because the accelerate d particles inside of the gap contribute to the 
total radiation emission outside the outer gap (IWada fc Shibatal 120071 ). 



2.3. Light curve 



Since the general feature of the light curve such as the standing phase of the pulse is 
mainly affected by the geometry of the emission regions, we produce a theoretical light curve 
by assuming constant emissivity. Figure [5] shows the theoretical light curve for a pulsar with 
the inclination angle a = 50°, at = 0.97, viewing angle Co — 76° and the azimuthal extension 
of the outer gap A<p p = 250°. The pitch angle at Rl is treated as a fitting parameter and 
is assumed to be sm<p(R L ) = 0.04. These parameters are chosen so that the modeled light 
curve explains the general features of the observation such as two peaks in a single period 
with a phase separation of ~ 140° between the two peaks. The breakdown of the light 
curve to show the contribution from the two poles separately is given in Figure El The color 
scheme is the same as the one for the emission pattern in Figure [2j As we shall see later, 
the inclination angle and the viewing angle chosen to explain the observed light curve also 
produces phase-resolved spectra which are consistent with observation. 



3. Energy spectra of the observed photons 
3.1. Acceleration and emission in the gap 



We adopt the local electric field equation in the CHR model for the region beyond the 
null charge surface. By assuming that the local electric field is decreasing in a quadratic form 
for the region between the null charge surface r nu n and the inner boundary of pair production 
region r in , we have 

r > r nu n, 



nB(r)h 2 (r) 



E» 



cs(r) 



r < r nuU . 



(7) 



For r > r nu n, E\\ is the vacuum solution given in ICheng. Ho fc Rudermanl (1l986al ). The 
vacuum solution is a good approximation for pulsars with thin gaps like the Crab pulsar. 



The radial distance r to the null charge surface varies with the field lines and the local 
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curvature radius also depends on the field lines. Therefore, strictly speaking, the electric 
field component along a magnetic field is a function of both the radial distance r and the 
polar cap angle <fr p , that is, En = E\\(r,(f) p ). 

When a relativistic charged particle which is accelerated continuously along a magnetic 
field line by the strong local electric field radiates by means of curvature radiation, the power 
gained by the accelerated particle as it goes through the electric potential, eE\\(r, </> p )c, is 
transformed to the power radiated as curvature radiation in order to maintain the equilib- 
rium. The total radiated power for each particle is l cur {r) = 2e 2 cy^(r)/3s 2 (r). The local 
Lorentz factor 7 e of the primary particles can be found by requiring eE\\(r, (p p )c = l cur , and 
hence, 

le{r) 



3s 2 (r) „ , \ 
__ii S|| ( rj ^) 



The characteristic energy of the radiated curvature photons is given by E cur (r) = f^-c^y- 
and Figure [7] shows the variation of E cur (r) along the field lines with the polar cap angle 
of P = 0° (solid-line), 90° (dashed-line), 180° (dotted-line) and 270° (dashed-dotted line). 
From Figure [7J we find that the particles are accelerated up to the ultra-relativistic regime 
so that 10 GeV photons are emitted in the outer gap accelerator of the Crab pulsar via 
curvature process. 



3.2. Synchrotron radiation and inverse Compton scattering from secondary 

pairs 

Since the relation E cur (r)dN 7 /dt ~ l cur {r)N is satisfied, the radiation spectrum of the 
primary particles in a unit volume is approximately described by 

d ( d 2 N 1 \ l cur {r)n 1 . . 

* 1 — (9) 



dV \dErydt J E CUT {r) E 7 

where n = dN/dV is the number density of the primary particles. Here, we use n = 
QB(r)/27rec which comes from the local Goldreich- Julian number density disregarding the 
angle between the local magnetic field direction and the rotational axis. Since the energy 
of each photon comes from curvature radiation, E~ < E cur (r). These primary curvature 
photons collide with the soft photons produced by synchrotron radiation of the secondary e + 
or e~ to produce even more secondary e 1 * 1 pairs. In this way, the synchrotron photons become 
abundant and nearly all the curvature photons are converted into secondary e 1 * 1 pairs. As a 
result, the energy distribution of the secondary e + or e~ is 

dN(r) „ 1 ( Em ^ r) d*N 7 (E 7 = 2E' e ) M 



'IE, " E,.h- dE.dt " ^ dEe 
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1 l cur (r)nB(r) 
E P 2iTecE cur (r) 



AV(r) In 



E,. 



E„ 



for t. 



77 



OO. 



(10) 



where the upper integration limit E max {r) is taken to be E CUT {r)/2 as the energy of each 
photon is divided between an e 1 * 1 pair, E e represents the energy of each e + or e~ and E e = 
2e A B 2 {r) sin 2 P(r)E 2 /3m^c 7 is the synchrotron energy loss rate. r 77 (i? 7 , r) is the attenuation 
depth for the absorption of the curvature photons. It will be very large in most of the 
magnetospheric regions except in the phases Trailing Wing 1, Bridge and Leading Wing 2. 
Detailed discussion on the absorption of curvature photons will be given in Section 

In the computer program, we divide the polar cap into N B equal divisions so that each 
division spans an angle of A<p p = 360°/Nb- Therefore, each division, represented by a 
magnetic field line, will take up a magnetic flux of ^Gap/Ns- The volume element AV{r) 
is represented by AA(r)Al(r) where AA(r) and Al(r) mean the area and the length of the 
tube-like volume element along a field line where observable photons are produced. Due to 
magnetic flux conservation, B{r)AA{r) = §go,p/Nb such that Equation fflOl) becomes 



dN(r) 
~dE~7 



1 l cur (r)Q® Gap Al(r) ^ ( E 

< 



E e 



E e 2necE cur (r)N B 
where $ Gap = 2nf(R L )B(R L )R 2 L . 

The photon spectrum of the synchrotron radiation is given by 



(11) 



E S yn(E^, r) 



dN(r) 
~dE~7 



N B 

E 

i=i 

V / 3e 3 i?(r) sin ip(r) 
mr,c 2 h 



dE^dt 



dE P 



syn 



E- 



N B 

E 

1=1 



dN{r) 
dE, 



E{x)dE e 



(12) 



where F(x) = x Ks (£)d£, Kb is the modified Bessel function of order 5/3, x — E 1 j E, 



syn \ 



and 



E S yn (^) 



3^5 



3heB(r) sin ip(r)E 2 /2mic 



(13) 



is the critical synchrotron photon energy. For the integration, the upper integration limit 
E ma x = E max (r) = E r/ur (r)/2, while the lower integration limit is chosen in a way that 
E min /m e c 2 = 20 iTakata. Chang fc Chenell2007h . 

Similarly, the photon spectrum due to inverse Compton scattering is given by 



(E„r) 



N B 

E 

i=l 



dN(r) 
dE P 



d 2 N 1 



dE e . 



(14) 
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If EpJrripC 2 » 1, the spectrum of t he inverse Compton scattered photons per electron 
(cf. §2.7 in iBlumenthal fc Gouldl . Il97ph is 



d 2 N 1 
dEydt 



ICS 



n 



syn 



[e,r) + n x (e,r)) 



x 



A(E e /m e c 2 ) 2 e 
2glng + (l + 2g)(l-g) + 



de 



(rg) 2 (l - q) 
2(1 + Tq) 



(15) 



where q = E 1 /T(l — Ei), T = 4(e/m e c 2 ) [E e j m e c 2 ) and E x = E 1 / E e . There are two possible 
soft photon sources for the inverse Compton scattering. The thermal photons with typical 
temperature T produced by the stellar surface produce a dominant absorption effect for 
photons emitted from locations close to the star. The number density of this kind of photons 
is given by 



Ux ^ r) = n^hcreMe/kT)-l {v) " (16) 
The other source of soft photons arises from synchrotron radiation. The absorption due to 
this kind of photons is more significant near the light cylinder. The number density of the 
synchrotron photons, n syn (e, r) is described by 



Hsyn f) 



syn 



cr 2 A£l(r) 

where AQ(r) is the solid angle of the beam of synchrotron photons and is estimated as 



(17) 



AQ{r) 



2t 



tp(r) 



i sin 9d6 



7np 2 (r) 



E„ in 



dN(r) 



in Eq. ([Tip is the sum of both the synchrotron energy loss rate and the inverse 

a T (m e ckT) 2 



Compton energy loss rat e due to thermal photons giv en by E e = ^^3" ' (hi(47 e fcT/ m e c 2 ) — 
5/6 — 0.5772 — 0.5700) (IBlumenthal fc Gouldl 119701) . The surface temperature of the Crab 
pulsar is taken to be 2 x 10 6 K (jTennant et al. 12001 ). 



In order that Equation f[T5|) is valid, q must be positive, so E\ < 1. Moreover, since the 
number of photons cannot be negative, we require the value in the bracket of Equation f|T5l) 
to be greater than or equal to zero. The upper integration limit €2 is chosen in such a way 
that E syn (e2, r) is very small for that particular field line; whereas the lower integration limit 
e\ is taken to be > 1 eV. 



3.3. Absorption of curvature photons 



In our previous studies (Cheng, Ruderman & Zhang 2000), we have anticipated that 
most of the curvature photons will be converted into the secondary pairs through the pair- 
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creation process with the magnetospheric X-rays. The typical pair-creation mean free path, 
which is estimated from l~ 1 (Ri c /2) ~ (1 — cos 9 co i)n x (Ri c /2)a^, where n x (Ri c /2) is the 
typical non-thermal X-ray number density, er 77 ~ cr^/3 is the pair-creation cross section 
with a T being Thomson cross section, and 6 C0 i(Ri c /2) ~ ^j2f(R lc /2)Ri c / s(Ri c /2) ~ \/0?2 
is the typical collision angle between the magnetospheric X-ray and the 7-rays emitted 
in the gap. For the Crab pulsar, the typical number density of X-rays is n x ~ L x (< 
E x >)/5Q(Ric/2) 2 c < E x >~ 6 x 10 17 cm 3 , where we used the typical energy < E x >~ 
(2m e c 2 ) 2 /(l — cos# co ;)/10 GeV ~ 260 eV, the typical non-thermal X-ray luminosity L x ~ 
5 x 10 34 erg/s, and the solid angle 5fl = 1 radian. As a result, the mean free path becomes 
/ ~ 2 x 10 7 cm (~ Ric/5) for the 10 GeV photons, and therefore we have believed that most 
of the 10 GeV photons emitted are converted into the secondary pairs. In the present paper, 
on the other hand, we take into the fact that some curvature photons emitted far away from 
the star (or near the light cylinder) could avoid the photon-photon pair creation process due 
to shorter escape distance and lower photon density near the light cylinder. In the present 
paper, we explicily calculate the optical depth for photons emitted at a given position. We 
can show that indeed some 10 GeV photons emitted near the light cylinder and almost all 
photons with energy below GeV can escape from the magnetosphere because their mean-free 
path becomes longer than the light radius. These survival curvature photons will contribute 
to the high energ y peaks around 1 GeV in the Trailing Wing 1 , Bridge and Leading Wing 2 
phases. Following iDing fc Cheng! (119971 ). the photon spectrum of the surviving 7-ray photons 
is obtained by 



TP — TP p -T{E^,r) 

- 1 cur,sur - 1 cw° 



(19) 



where F cur is the curvature spectrum and is described by 

F (E r) = dN ^' F(x) 
cur[ 7 ' ] 2nhsE^ 1 J 



(20) 



where dN is the number of primary e 1 * 1 pairs in the emission region and is give n by dN 



n^jAAAl and x = E y / E cur {r). The attenuation depth t(E 7 , r) is calculated from (IJauch fc Rohrlich 
1976h 



r{E 7 ,r) =l(r) 



[n syn (e, r) + n x (e, r)]a 77 (E 7 , e)de 



(21) 



with Z(r) being the distance between the emission location and the light cylinder and 
<J 77 (i? 7 ,e) being the cross section for photon-photon pair creation and is given by 



<r 77 (£ 7 ,£) = — a T {\ 
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(3 - v A ) In 



1 + v 



2v(2-v 2 ) 



(22) 
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where 

/ m 

(23) 



E^E 



The total photon flux received on Earth is 

F{EJ = Yl i F syn(E 7 , r) + F ICS (E^ r) + F c „ r (E„ r)] (24) 

r 

where D = 2 kpc is the distance of the Crab pulsar from the Earth and AQ is the solid 
angle chosen to be 1 sr for simplicity. 



3.4. Phase-resolved spectra 

The computation of the photon flux is divided into four steps. Firstly, we calculate the 
synchrotron radiation photon flux. Then, we use this synchrotron flux to calculate the inverse 
Compton scattering photon flux. Next, we adjust the peak intensities of the synchrotron and 
the inverse Compton scattering spectra by requiring J E 1 F{E 1 )dE 1 for synchrotron spectrum 
alone before the consideration of the inverse Compton scattering radiation to be the same 
as the sum of the synchrotron and the inverse Compton scattering spectra. However, the 
relative peak intensities of the synchrotron and the inverse Compton scattering spectra are 
kept unchanged. We employ this third step to comply with energy conservation because 
all the energy in inverse Compton scattered photons comes from the synchrotron photons. 
Finally, we calculate the survival curvature radiation. 

Figure [10] shows the observed data of the phase-averaged spectrum of the Crab pulsar 
and the theoretical fitting spectrum with f(Ri) = 0.2 and sin ^{Rl) = 0.06, calculated 
by using the synchrotron self-Compton mechanism together with the survival curvature 
photons from 100 eV to 10 GeV. The blue and red lines show the emission beyond the null 
charge surface (i.e. from the south pole) and between the inner boundary and the null charge 
surface (i.e. from the north pole), respectively. For each colour, the dashed line represents the 
synchrotron spectrum, the dotted line represents the inverse Compton scattering spectrum, 
the dash-dotted line represents the survival curvature spectrum and the black solid line is 
the sum of the three spectra from both poles. The fitting spectrum explains the observation 
from 100 eV to 10 GeV. 



Figure [TT] shows the predicted phase-resolved spectra as a break-down of the phase- 
averaged spectrum in Figure [10J The colour and line sch emes are the same as in Figure [TOl 
The phase intervals are defined in the same way as in iFierro et al.l (119981 ) . The pulse of 
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the Crab pulsar is divided into 8 phases, namely Leading Wing 1, Peak 1, Trailing Wing 
1, Bridge, Leading Wing 2, Peak 2, Trailing Wing 2 and Off Pulse. The criteria for the 
division is listed in Table [TJ We choose the azimuthal angle of 52° from the Q-/J, plane to 
be phase zero. The footprints corresponding to each phase are plotted for both poles in 
Figure [H Although our prediction is not too good (for example, in LW1, the peak at about 
100 MeV cannot be produced; whereas in P2 and TW2, the peaks at about 100 MeV are 
much too high), most of the spectral features can be explained. However, it has been argued 
that the north pole and the south pole need not to be perfectly sym metric. Since the dipole 
may not be at the centre of the star, just like the sunspot geometry (jRudermanl Il99ll ) . it is 
possible that the two poles could have small differences. Furthermore, sin <p(Rl) is calculated 
according to a simple static dipolar field, which is clearly not a good approximation in the 
outer magnetosphere of a realistic rotating dipolar field with current flow. In the best fit, 
we thus allow the gap sizes to be different for the two poles and the pitch angles to vary for 
individual phases. 



Figure [T2l shows the comparison between the observed data of the phase-resolved spectra 
of the Crab pulsar and the theoretical best fitting results. Table [2] summarizes the values of 
/(Rl) and sm(p(R L ) for different phases for each pole. We find that the fitting pitch angles 
for the phases consisting of field lines at a polar angle closer to (f) p = 0° are smaller, while 
those for the phases consisting of field lines at a polar angle closer to <ft p = 180° are larger. 
This agrees with the theory that the outer gap is thinner around <p p = 0° and thicker around 
(p p = 180°. 

As we consider the emission region to extend inwards to an inner b oundary inside the null 
charge surface, we can fit the spectrum of LW1. This cannot be done in lCheng. Ruderman fc Zhang 
( 120001 ) which considered the emission region to be from null charge surface to light cylinder 
alone, i.e. one pole only. Figure [T21 shows that this phase is mainly contributed by the radia- 
tion fr om the n orth pole. Another phase that canno t be produced in lCheng. Ruderman fc Zhang 
tood ) is TW2. ICheng. Ruderman fc Zhan^ (boooh estimated A0 P ~ 160°, but in Section^ 
we estimated that A(j) p can be ~ 250° for the extended new emission geometry. The widening 
of the azimuthal extension of the outer gap allows us to obtain a reasonable fitting for the 
phase TW2. 



From Figure [12j we can see that except LW1 and the 100 MeV regime of PI, all the 
phase-resolved spectra are dominated by the radiation from the south pole. In general, the 
spectra are mai nly contributed by the part from the null charge surface to the light cylinder, 
so the results in ICheng. Ruderman fc Zhang] (120001 ) are good approximations. The peaks at 
the high energy regime near several GeV in TW1, Bridge and LW2 are believed to be the 
survival curvature radiation emitted at locations far away from the star (cf. Figure [9]). This 
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agrees with the model that most of the curvature radiation has been converted into e pairs 
and then emitted as synchrotron and inverse Compton scattering photons. 

For the 3 phases between the peaks, i.e. in the middle row of Figure [T21 namely TW1, 
Bridge and LW2, the synchrotron peaks from the south pole are much wider than their 
counterparts from the north pole. In fact, they show a kind of double-peak feature with a 
plateau around 1 MeV. According to Figure the radiation of each of these 3 phases are 
actually emitted from two well-separated regions. For TW1, the radiation is emitted from an 
outer region around l.lbRi < r < 1.7Rl and an inner region around 0.46-Rl < r < 0.67Rl- 
For Bridge, they are 1.28R L < r < 1.72R L and 0.31R L < r < 0A5R L . For LW2, they are 
1.15.Rl < r < \.TJR L and 0.29-R^ < r < 0.31Rl- In order to know more clearly this feature, 
we plot a further decomposed phase- resolved spectra for TW1, Bridge and LW2 in Figure [TBI 
In Figure [13], we decompose each blue line in Figure [12] into two components. The green 
lines represent the radiation coming from the outer emission region and the magenta lines 
represent those from the inner region. 

From Figure [13] we can notice that the radiation emitted at a larger r, i.e. from the 
outer region, will produce a synchrotron peak at a lower energy. This can be explained by the 
typical synchrotron photon energy given in Equation (fT3l) with E syn oc B(r) simp(r) oc r~ 2,5 . 
Let us take the Bridge as an example. We choose the representative r to be the midpoint of 
the emission region from Figure [9] i.e. 1.25Rl for the outer region and 0.35Rl for the inner 
region. In this case, the ratio of r is about 3.57 and hence the characteristic E syn of the 
inner region is about 24 times that of the outer region. From the modelled fit, the ratio of 
the peaks is about 27. Therefore, the theory matches quite well with the observation in this 
case. 



Conclusion and discussion 



We have used the modified three dimensional outer magnetosphere gap model with 
the inner boundary of the outer gap being extended from the null charge surface to near 
the stellar surface. The exact location of the inner boundary does not affect the fitting 
results. The "inwardly-extended" part of the outer gap contributes to LW1 and TW2 of 
the light curve with a slight modification of PI and P2. Such modified outer gap geometry 
also plays a vital role in expla ining the optical polarizatio n properties of the Crab pulsar 



flTakata. Chang fc Chendl2007f ) . Together with the results of iTakata. Chang fc Chend (120071 ) 
on the gap accelerator, this model explains the pulse profile, phase-resolved spectra and 
polarization of the Crab pulsar. Also, the outer gap model can expla in the observed comple x 
morphology change of light curve as a function of the photon energy (ITakata fc Chang||2007l ). 
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Four adjustable parameters are used to simulate the light curve: the inclination angle 
of the magnetic axis to the rotational axis a, a 1; the viewing angle also to the rotational axis 
( and the emission cone pitch angle due to geometry of the emission location sin <p(Rl). As 
constrained by the phase separati on of th e doub le peaks, we choose the values a = 50° and 
( = 76°. From radio observations. iRankinl (119931 ) estimated a ~ 84° with ( unknown. By us- 
i ng th e polarimetric observations at frequencies between 1.4 and 8.4 GHz. iMoffett fc Hankins 
(119991 ) calculated a ~ 56° and ( = 117°. Therefore, further observations is required in order 
to determine these two parameters more accurately. 

For the photon spectra, our model fitting requires J(Rl) ~ 0.2, which is larger than the 
theoret ical estimation ( f(Rr,) ~ 0-12) by a factor of ~ 2. However, the theory (jZhang fc Cheng 



19971 ; ICheng. Ruderman fc Zhang||2000l ) has assumed a vacuum dipole potential. It has been 



pointed out that in order to explain the observed radiation power on high energy region, 



Hirotani 


2006a. 


2007; 


Takata et al. 


2006) 



should be larger than the vacuum one. 

Moreover, we have taken the stellar radius of a neutron star to be 10 6 cm. Due to the 



lack of a direct method to determine the size of a neutron star ( ILattimer fc Prakashl 120041 ) . 
and the equation of state inside a neutron star given by different theoretical models do not 
provide a unanimous value for the neutron star size, we can only determine the magnetic 
moment, i.e. B p R? , of the pulsar from the energy loss rate. However, although the dipole 
radiation comes from the spin down power of the star, it is not all of it. For example, the 
vacuum formula of the dipole radiation for the aligned rotator is usually used t o infer the 
dipole moment of the pulsar. On the other hand, MHD solution (jSpitkovskyll2006l ) suggested 
that the effect of the current increases the spindown luminosity from the the vacuum formula. 
Therefore, we may have overestimated the value of the magnetic moment and therefore the 
strength of the magnetic field. Since E syn oc B(r) hence B p R 3 , the synchrotron peak may 
have been overestimated. Similarly, the peak synchrotron flux F syn oc B(r) may also be 
overestimated. Since the power of synchrotron radiation and inverse Compton scattering 
can be compared by the ratio of the local magnetic energy density and the photon energy 



density, 



Psyn ^ Ub ^ B\r)l2-K 

PlCS U-y E s yn(r)n s yn(Esyn,r) 



which is independent of B P R? shows that Fjcs may 



not be affected. On the other hand, since F syn may have been overestimated, n syn should be 
smaller. Hence less curvature photons should be absorbed. 

The Crab pulsar is one of the most important pulsars because of its close distance and 
strong radiation signals in all frequencies. Although observation data of it have been gath- 
ered for about 40 years, the underlying physics involved in this pulsar are not completely 
understood. For example, the structure of the charged particle accelerator in the magneto- 
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sphere can be d ifferent in the polar c ap model, the outer gap model and the newly proposed 



caustic model (iDyks & Rudakl 120031 ) which proposed a two pole model with the thin gap 
regions starting from each polar cap to the light cylinder. We hope that further observation 
by the more sensitive satellites in the future missions, e.g. GLAST, can bring more crucial 
and critical information so that we can determine the mechanisms that are truly taking 
place. For example, as shown in Figure [T2], we predict that there is a clear component (30 
MeV - 30 GeV), which is mainly dominated by the survival curvature photons, in TW1, 
Bridge and LW2. Our model also suggests that the high energy turn over in TW1, Bridge 
and LW2 is ~ 10 GeV, whereas the high energy turn over of LW1 is only ~ 100 MeV (much 
lower than the other phases). This difference results from the consequence of the radiation 
in LW1 being emitted from a different pole in the region between the inner boundary and 
the null charge surface, where the local electric field is much weaker (cf. Eq. [7]). With an 
increase in sensitivity and a widening of the energy range than EGRET, GLAST may be 
able to confirm our predictions. 
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Fig. 1. — The 3-dimensional view of the last closed field lines for a = 50° in the rotating 
frame. The rotational axis is the z-axis. The f2-/x plane is the x-z plane. Crab period is 
used: P = 0.33 s, so R L ~ 159i? w 1.59 x 10 s cm. 




Fig. 2. — Photon emission pattern for a = 50° and a\ = 0.97. 
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Fig. 3. — Close-up illustration of two field lines that approximate two concentric circles. 
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Fig. 4. — The radial distance to the inner boundary and the null charge surface on the last 
closed field lines. 4> p refers to the polar cap angle. The angle P = 0° represents the Q-fi 
plane. 
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Fig. 5. — Theoretical light curve for the Crab pulsar. The fitting parameters are a = 50°, 
ai = 0.97, Co = 76° and smip(R L ) = 0.04. 
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Fig. 6. — The breakdown of the light curve in Figure [5] to show the contribution from each 
pole. 
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Fig. 7. — Curvature photon energy vs r for the field lines at P = 0°, 90°, 180° and 270° 
represented by solid, dashed, dotted and dash-dotted lines, respectively. 
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Fig. 8. — Footprints that correspond to particular phases. The upper panel is the north pole 
and the lower panel is the south pole, a = 50°, ai = 0.97 and q = 76°. Some footprints do 
not belong to any phase and some of them belong to more than one phases. From (f> p = 73.5° 
to 223.5°, there is no radiation produced along the line of sight. 
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Fig. 9. — The emission locations for photons which have been counted in Figure El 
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Fig. 10. — Phase-av eraged spectrum of the Crab pulsar. The observed data are taken from 
Kuiper et all toQlh . The fitting parameters are a = 50°, a x = 0.97, q = 76°, f(R L ) = 0.2 
and sintp(RL) = 0.06. 
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Fig. 13. — The breakdown of the south pole spectra (blue) of the phases TW1, Bridge and 
LW2 in Figure [12j each being decomposed into two components according to the emission 
locations. 
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Table 1: Phase component definitions for the Crab pulsar. 



Component 


Abbreviation 


Phase interval 


Width 


Leading Wing 1 


LW1 


8.8° - 


30.4° 


21.6° 


(0.06) 


Peak 1 


PI 


30.4° - 


- 66.4° 


36° 


(0.10) 


Trailing Wing 1 


TW1 


66.4° - 


102.4° 


36° 


(0.10) 


Bridge 


Bridge 


102.4° 


- 142° 


39.6° 


(0.11) 


Leading Wing 2 


LW2 


142° - 


167.2° 


25.2° 


(0.07) 


Peak 2 


P2 


167.2° - 


- 206.8° 


39.6° 


(0.11) 


Trailing Wing 2 


TW2 


206.8° - 


- 239.2° 


32.4° 


(0.09) 


Off Pulse 


OP 


239.2° 


-8.8° 


129.6° 


(0.36) 
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Table 2: Summary of fitting parameters sm (p(Ri,). We have taken }\Rl) to be 0.25 for the 
north pole and 0.22 for the south pole. 

Phases sin <p(Rl) 

North Pole South Pole 



LW1 


0.02 




PI 


0.03 


0.08 


TW1 




0.08 


Bridge 




0.08 


LW2 




0.06 


P2 


0.08 


0.08 


TW2 


0.07 


0.08 



